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ABSTRACT 

We present the measured projected obliquity - the sky-projected angle between the stellar spin 
axis and orbital angular momentum - of the inner planet of the HAT-P-17 multi-planet system. We 
measure the sky-projected obliquity of the star to be A = 19l}g degrees by modeling the Rossiter- 
McLaughlin (RM) effect in Keck/HIRES radial velocities (RVs). The anomalous RV time series 
shows an asymmetry relative to the midtransit time, ordinarily suggesting a nonzero obliquity - 
but in this case at least part of the asymmetry may be due to the convective blueshift, increasing 
the uncertainty in the determination of A. We employ the semi-analytical approach of Hirano et al. 
(2011) that includes the effects of macroturbulence, instrumental broadening, and convective blueshift 
to accurately model the anomaly in the net RV caused by the planet eclipsing part of the rotating 
star. Obliquity measurements are an important tool for testing theories of planet formation and 
migration. To date, the measured obliquities of ~50 Jovian planets span the full range, from prograde 
to retrograde, with planets orbiting cool stars preferentially showing alignment of stellar spins and 
planetary orbits. Our results are consistent with this pattern emerging from tidal interactions in the 
convective envelopes of cool stars and close-in planets. In addition, our 1.8 years of new RVs for this 
system show that the orbit of the outer planet is more poorly constrained than previously thought, 
with an orbital period now in the range of 10-36 years. 



1. INTRODUCTION 

HAT-P-17 is an early K dwarf star that hosts a tran- 
siting Saturn-mass planet (planet b) on a 10.3 day orbit 
and a more massive outer companion (planet c) on a long- 
period orbit (Howard et al. 2012, hereafter H12). Tran- 
sits of planet b were discovered in 2010 by the Hungarian- 
made Automated Telescope Network (HATNet, Bakos 
et al. 2004). Followup Keck/HIRES RVs were used to 
measure the mass of planet b and enabled the discov- 
ery of planet c. More than 150 hot Jupiters have been 
discovered, but it appears that hot Jupiters tend to lack 
additional short-period giant planet companions (Stef- 
fen et al. 2012). HAT-P-17 is one of only six of sys- 
tems with a transiting Jovian-sized planet and an addi- 
tional substellar companion. The five other systems in- 
clude HAT-P-13 (Bakos et al. 2009), HAT-P-31 (Kipping 
et al. 2011), Kepler-9 (Holman et al. 2010), Kepler-30 
(Fabrycky et al. 2012), and KOL94 (Hirano et al. 2012). 
These rare multi-planet transiting Jovian systems pro- 



^ Institute for Astronomy, University of Hawaii at Manoa, 
2680 Woodlawn Dr, Honolulu, HI 96822 

^ Department of Physics and Kavli Institute for Astrophysics 
and Space Research, Massachusetts Institute of Technology, 
Cambridge, MA 02139, USA 

^ California Institute of Technology, Pasadena, CA, USA 
University of California, Berkeley, Berkeley, CA, USA 

^ Princeton University, Department of Astrophysical Sciences, 
Princeton, NJ, USA 

® University of Notre Dame, Department of Physics, Notre 
Dame, IN, USA 

^ Alfred P. Sloan Fellow 

* Packard Fellow 

^ Department of Astronomy and Astrophysics, 525 Davey Lab- 
oratory, The Pennsylvania State University, University Park, PA 
16802, USA 



vide important insight into the formation and evolution 
of hot Jupiters. 

Current theory suggests that Jovian planets form at or- 
bital distances of >1 AU where additional protoplanetary 
solids (ice) augment their formation. They then migrate 
inwards to become hot Jupiters. Popular theories that 
attempt to explain their resulting close-in orbits involve 
a 3rd body (in addition to the Jovian planet and it's host 
star) that perturbs the orbit of the soon-to-be hot Jupiter 
and excites high eccentricities through either the Kozai 
mechanism or planet-planet scattering. This highly ec- 
centric orbit then decays through tidal interactions into 
a close-in circular orbit (Nagasawa et al. 2008; Fabrycky 
& Tremaine 2007; Naoz et al. 2011). This scenario would 
produce hot Jupiters with a large range of orbital obliq- 
uities. Others suggest that hot Jupiters migrate within 
the circumstellar disk from which they formed through 
interactions with the disk (Lin et al. 1996). In this case 
we expect that all of the bodies would lie in coplanar 
orbits that are all well-aligned with the stellar spin axis. 
If the orbit of planet b is aligned to host star's spin, it 
would suggest that this system was formed by migration 
rather than perturbation if the two planets are copla- 
nar. A coplanar and apsidally locked geometry would 
also allow for a precise measurement of the interior den- 
sity structure of planet b (Batygin et al. 2009; Mardling 
2010). 

An emerging trend suggests that hot Jupiters around 
cool stars (Toff < 6250 K) with large convective envelopes 
tend to be better aligned with their host star's rotation 
axis (Albrecht et al. 2012). Tidal energy is most effi- 
ciently dissipated by turbulent eddies in the convective 
regions of stars (Zahn 2008). As a result, the rate of 
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tidal dissipation depends on the mass of the convective 
envelope. Strong tidal interactions with the convective 
envelope force the system into alignment in a relatively 
short time. Stars hotter than 6250 K have small or no 
convective envelopes, and it takes much longer for the 
system to align (Winn et al. 2010b). HAT-P-17 is a cool 
star with Toff ^ 5200 K, but planet b's orbital distance 
is relatively large making tidal interactions weak. Ac- 
cording to the tidal figure of merit devised by Albrecht 
et al. (2012), we would expect the tidal dissipation rate 
for this system to be too slow to cause obliquity damp- 
ing, despite the star's thick convective envelope. This 
makes HAT-P-17 an interesting test case. 

In this work we revisit the orbital parameters of planet 
c with new Keck/HIRES RV and Keck/NIRC2 adaptive 
optics images, and present a measurement of the sky- 
projected orbital obliquity of the star relative to planet 
b. In §2 we discuss our observational techniques. We 
discuss our RV and RM modeling and results in §3, and 
in §4 we interpret and summarize our findings. 

2. OBSERVATIONS 

2.1. Keck/HIRES Spectroscopy 

Since the publication of Howard et al. (2012), we have 
measured the RV of HAT-P-17 (V=10.54) for an addi- 
tional 1.8 years using HIRES (Vogt et al. 1994) on the 
Keck I telescope. Wc adopted the same observing strat- 
egy and Doppler analysis techniques described in Section 
2.3 of H12. In brief, we observed HAT-P-17 through a 
cell of gaseous iodine and measured the subtle Doppler 
shifts of the stellar lines with respect to the reference io- 
dine lines using a forward modeling analysis (Butler et al. 
1996). 

Our observations were designed to measure the Kep- 
lerian orbits of HAT-P-17b and c and also to measure 
the obliquity of HAT-P-17. For the latter, we observed 
a transit of HAT-P-17b on UT 26 August 2012. Our ob- 
serving sequence lasted nearly six hours and bracketed 
the 3.2 hour long transit. We made 42 observations of 
~500 second duration separated by 45 second detector 
reads. To constrain the Keplerian slope, we made three 
additional observations on the same night approximately 
3.8 hours after transit egress. 

Julian dates of the photon-weighted exposure mid- 
times were recorded during the observations, and then 
later converted to BJDtdb using the tools described in 
Eastman et al. (2010)^°. The photon- weighted exposure 
times are only accurate to ~1 second due to internal lim- 
itations of the exposure meter. 

The complete set of RV measurements and their un- 
certainties are listed in Table 1. These 100 RVs include 
42 RVs from Howard et al. (2012), 45 new RVs taken 
on the night of UT 26 August 2012 to measure the RM 
effect, and 13 additional RVs taken sporadically between 
2010 April and 2012 December to measure the orbit of 
HAT-P-17C. 

2.2. KECK/NIRC2 adaptive optics imaging 

In order to search for additional companions and 
sources of possible photometric dilution, we obtained 

IDL tools for time systems conversion; 
http://astroutils.astronomy.ohio-state.edu/time/ 



TABLE 1 
Radial velocity data'' 



Time RV (Trv 

BJDtdb - 2440000 m s^^ m s^^ 



14396.8272772 


-5.25 


1.62 


14397.7946382 


-32.79 


1.60 


14427.7815123 


-10.63 


1.58 


14429.8199962 


-68.92 


1.78 


14430.8485952 


-97.84 


1.90 


14454.7162454 


9.70 


2.66 


14455.7078194 


14.89 


1.90 



This tabic is available in its entirety 
in macliiuc-rcadablc form in the online 
journal. A jjortion is shown here for 
guidance regarding its form and content. 



high spatial resolution images of HAT-P-17 using NIRC2 
(instrument PI: Keith Matthews) at the Keck II tele- 
scope on 2012-05-07 UT. Photometric dilution would af- 
fect the radius of HAT-P-17b measured by H12, and the 
presence of a physically associated companion would put 
constraints on our radial velocity fit. In addition, a sta- 
tistical sample of the wide companions to exoplanet host 
stars may help our understanding of planetary forma- 
tion mechanisms. Our observations consist of dithered 
images acquired using the K filter (central wavelength 
= 2.12/xm). We used the small camera setting to pro- 
vide fine spatial sampling of the instrument point spread 
function. The total on-source integration time was 16.2 
seconds. Images were processed using standard tech- 
niques to flat-field the array, replace hot pixels, subtract 
the thermal background, and align and co-add individual 
frames. 

3. ANALYSIS 
3.1. Radial velocities 

With 1.8 years of new radial velocities we revisited 
the orbital parameters of the outer companion in the 
HAT-P-17 system (planet c). We analyzed the entire 
dataset with a custom version of EXOFAST^^ (East- 
man et al. 2012) ported to Python (ExoPy hereafter). 
ExoPy utilizes the Differential-Evolution Markov Chain 
Monte Carlo (DE-MCMC) technique (Ter Braak 2006) 
to find the best-fitting parameters and their associated 
uncertainties. We also ported a subset of the RVLIN^^ 
(Wright & Howard 2009) package to Python for quick 
calculation of the Keplerian orbit model. We fit a 12 pa- 
raniciter model to the radial velocity data that included 
the period {Pb), time of transit (Ttra,b), eccentricity (ef,), 
argument of periastron of the star's orbit (wb), radial 
velocity semi-amplitude {Kb) of planet b, the same pa- 
rameters for planet c (Pc, Ttra.c, "^c, Kc)^ the center 
of mass velocity of the system normalized to an arbitrary 
reference spectrum (7), and a stellar "jitter" term. 

We computed 24 DE-MCMC chains in parallel, con- 
tinuously checking for convergence following the pre- 
scription of Eastman et al. (2012). We considered the 
chains well-mixed and halted the DE-MCMC run when 
the number of independent draws (T^, as defined in 
Ford 2006) was greater than 1000 and the Gelman-Rubin 

IDL code available at; http://astroutils.astronomy.ohio- 
state.edu/exofast/ 

IDL code available at; http://exoplaiiets.org/code/ 
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Fig. 1. — Top: Keck/HIRES RV measurements for HAT-P-17 as a function of BJDxdb with the best-fitting two-planet model and 
associated residuals found by minimization of the DE-MCMC chains. A stellar "jitter" term (see Table 2) has been added in quadrature 
to the measurement errors. Data taken during the transit of HAT-P-17b for the purpose of measuring the RM effect was excluded from 
the RV fit and are not included in this plot. Data to the right of the vertical dashed line are new to this work, and data to the left are from 
H12. Bottom: same RV measurements phase-folded to the orbital ephemerides of planets b (upper) and c (lower). Phase corresponds to 
the time of mid-transit (or hypothetical transit). In each case the orbit of the other planet and an arbitrary center of mass velocity relative 
to a template spectrum (7) has been removed. 

The results of the RV analysis are presented in Table 
2. All of the parameters for planet b are consistent with 
the values from H12. However, with the new RV data we 
can now see that the period of planet c is much longer 
than the initially reported period of 1620 ± 20 days. 
Our RV time-series span a total timespan of 1869 days. 
We still have not seen a complete orbit of planet c, and 
thus the fit is quite poorly constrained. Figure 2 shows 
the probability distributions of Mc sin ic and eccentricity, 
and indicates that the allowed mass range (Su) for planet 
c is 2-6 Mj. 

We also explored the possibility of a 4th body causing 
a linear trend (7) in the radial velocities in addition to 
the signal from planet c. Fits that included 7 as a free 
parameter preferred a slope consistent with zero (7 — 
— S.T^s'o s~^yr~^). To assess the validity of adding 
one more free parameter to our model we calculate the 
bayesian information criterion (BIG): 



statistic (Gelman et al. 2003; Ford 2006; Holman et al. 
2006) was within 1% of unity for all parameters. In 
order to speed convergence, ensure that all parameter 
space was adequately explored, and minimize biases in 
parameters that physically must be finite and positive, 
we step in the modified and/or combinations of parame- 
ters shown in Table 2. Namely, due to the highly corre- 
lated uncertainties of e with uj and v sin i^^ with A we step 
in Y^cosw, "v/esinw, ^/v sini* cos A, and ^Jvsvaii, sin A 
(Eastman et al. 2012; Albrecht et al. 2012). 

We assigned Gaussian priors to Pf, and Ttra from the 
values given in the HAT-P-17b,c discovery paper {Pb = 
10.338523 ± 9 X 10"^ days, Ttra.b = 2454801.16943 ± 
2 X 10^"* BJDtdb, H12) that came from the highly con- 
straining photometric transit data, and we assigned uni- 
form priors to all other step parameters. We ignored any 
variation in transit times due to perturbations caused by 
planet c, but these are expected to be at least an order of 
magnitude smaller than the propagated uncertainty on 
Ttra,b (H12). We excluded the spectroscopic transit data 
from the Keplarian radial velocity fit because the higher 
density of data on the one night could bias the results 
in the presence of short-term (~hours) trends. The best 
fitting values and upper and lower "Icr" errors for each 
parameter were determined by taking the median, 85.1, 
and 15.9 percentile values, respectively, of the resulting 
posterior distributions. 



BIG = x +k\nn 



where k is the number of degrees of freedom, and n is 
the number of data points in the fit (Liddle 2007) . The 
BIG increased when 7 was a free parameter (104 vs. 100 
with or without 7 as a free parameter respectively). The 
BIG increase, our model fit favoring 7 = 0.0, and the 
AO image (see §3.4) all indicate that the data are better 
described by a model with 7 fixed at zero. 



4 



Fulton et al. 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 





TABLE 2 
Radial velocity MCMC results 




2 3 4 5 6 
Msin/[Mj] 



Fic;. 2. Probability distribution of eccentricity vs. minimum 
mass (M sin t) in Jupiter masses for planet c from the DE-MCMC 
analysis. The dashed lines are 68%, 95%, and 99% confidence 
intervals and the white dot is the median value. The median value 
is offset from the mode (most likely value) due to the asymmetric 
posterior distributions. 

3.2. Spectroscopic transit 

At first glance the RM data follow the typical red- 
shift then blueshift pattern of a spin-orbit-aligned sys- 
tem. However, the data do not cross zero until slightly 
after the time of mid-transit. The small asymmetry in 
the RM curve (Figure 3) suggests a slight misalignment. 
We also used ExoPy to analyze the spectroscopic transit 
data. Our model of the RM effect takes the form of 

iiMnet = Av + Vcb + 8(1 - Tt,a) + vcm (2) 

where Av is given by equation 16 of Hirano et al. (2011) 
and is discussed in more detail in section 3.2.1 below. 
Vcb is the anomalous radial velocity shift due to the 
convectivc blueshift (Shporer & Brown 2011, discussed 
in section 3.2.2 below), S is the radial velocity slope ob- 
served during transit due to the orbital motion of HAT- 
P-17b-|-c, t are the flux-weighted exposure mid-times of 
the observations in BJDtdb, T^ra is the BJDtdb of mid- 
transit, and Vcm is an arbitrary additive constant veloc- 
ity. Pb is constrained to the value obtained in the RV 
analysis, and Ttia.b is constrained by propagating the 
error on Ph and Ttra,b found from the Keplcrian anal- 
ysis to the transit epoch of the night of 2012 Aug 26 
(Ttra.b = 2456165.8553 ± 0.0012 BJDtdb)- The ampli- 
tude of the HJDuTC to BJDtdb correction applied to 
the RV data was ^67 seconds, or about half of the prop- 
agated uncertainty on the mid-transit time which high- 
lights the importance of working in a standardized and 
consistent time system. S and vcm were assigned uni- 
form priors in order to account for any trends present 
during the observed transit. 

3.2.1. Semi-analytical Rossiter-McLaughlin model 

The shape and amplitude of Av depends on nine pa- 
rameters. Five describe the decrease in flux as the planet 



Parameter 



Value 



Units 



1.0144585 ±3.7e - 07 
2454801.1702 ±0.0003 



0.5442 



R.V Step Parameters: 
log(P6) 

log(Pc) 
^tra,c 
•^/e^ COS U!c 

^/eZ sin u>c 

7 
7 
7 

logQitter) 



+0.0052 



0.0051 

-0.214 ±0.016 

1.7678 ±0.0051 
, +0.38 
'^■''^ -0.2 

2454146 tiro 
-0 63 +"1^ 

-0 017 +° °^* 

-0.061 

1 fiSQ +008 
l.bSy _o.o61 

20 +2'' 
-16 

= 0.0 

= 0.0 

n QIO +0.076 



log(days) 
BJDtdb 



m s ^ 
log(days) 
BJDtdb 



m s~^ 
m s~^ 
m s~^day~^ 
m s~^day~^ 
log(m s~^) 



RV Model Parameters: 
10.338523 



Pb 

^tra,b 

Sb 

l^b 

Kb 
Pc 

'^tra,c 

Be 

U>c 

Kc 
7 

7 

7 

jitter 



+ 8.8e-06 
-8.9e-06 

2454801.1702 ±0.0003 
0.3422 ±0.0046 
201.5 ±1.6 
58.58 

+7700 
-2100 
-100 
-170 
0.23 
-0.17 
+5.3 
6.7 
+9.9 
6.4 
+27 
-16 

= 0.0 

= 0.0 

^•UO _o.35 



5584 
2454146 
0.39 
181.5 
48.8 
20 



days 
BJDtdb 

degrees 
m s~^ 

days 
BJDtdb 

degrees 
1 



m s 
m s 

m s~^day 

m s^ 



m s 



1 
1 

-2 
1 



RV Derived Parameters: 




Mf, 


n t^Qo +0.018 

U.JOZ _g 017 


Mj 


a-b 


0.0882 ±0.0014 


AU 


Mc sin ic 




Mj 


Oc 


5.6 


AU 



transits its host star; the planet to star radius ratio 
{Rp/R^), the semi-major axis of the orbit in units of 
stellar radii {a/Ri,), the inclination of the orbit relative 
to our linc-of-sight (i), and two quadratic limb darkening 
coefficients (mi,U2). We assigned Gaussian priors to Pb, 
7tra,bj Rp/ R*y 0,/R*, and i from the values given in H12 
as these are poorly constrained by the RM data alone. 

Two more geometrical parameters contribute to the 
shape of the spectroscopic transit; the rotational velocity 
of the star projected onto the plane of the sky {v sini^,), 
and the angle between the rotational axis of the star 
projected onto the plane of the sky and the planet's or- 
bital angular momentum vector (A). We adopt a value 
of vsinii, = 0.3 ± 1.5 km s~^ as a Gaussian prior that 
was obtained from the Spectrocospy Made Easy (SME, 
Valcnti & Piskunov 1996) analysis performed in H12. 

Some of the orbital parameters from the RV analysis 
also have a small effect on the timing and duration of the 
transit. We assigned Gaussian priors to and ojb from 
the results of the RV analysis. 
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Fig. 3. — Left: RV variation during the transit of HAT-P-17b due to the RM effect with best-fit model found from x minimization 
overplotted. The upper panel includes the RV variation due to the orbital motion of HAT-P-17b, the middle panel shows the data and 
model with the orbital motion removed, and the bottom panel shows the residuals to the model. The solid blue line represents our adopted 
model including all line broadening effects and the convective blueshift, the dashed green line shows the model without the convective 
blueshift, and the dot-dashed magenta line shows an idealized model in which the line profiles are described only by rotational broadening, 
similar to the approach of Ohta et al. (2005). Note the asymmetry in the RM curve caused by a combination of the convective blueshift 
and a slight misalignment. Right: Posterior distribution of tisinu vs. A from the DE-MCMC analysis of the spectroscopic transit. The 
dashed lines are 68%, 95%, and 99% confidence intervals and the white dot is the median value of the distribution. The dot-dashed line in 
the right histogram shows the 0.3 it 1.5 km s~^prior on ii sin j from the SME analysis of H12. 

The convective blueshift {Vcb) is caused by the net 
convective motion of the stellar photosphere. Hotter 
material from below the photosphere rises upward to- 
wards the observer due to convection and is only par- 
tially canceled by downwelling cold material, causing a 
net blueshift of order 1 km s~^. Since we are only in- 
terested in relative radial velocities this net blueshift is 
unimportant. However, because the convective blueshift 
is strongest near the center of the star and weaker near 
the limbs, the transiting planet occults areas of the star 
that have different contributions to the net convective 
blueshift. This causes a time-varying component of the 
convective blueshift during the spectroscopic transit of 
order 2 m s^^ . Since the v sin of HAT-P-17 is relatively 
low, the amplitude of the spectroscopic transit signal is 
only about 7 m s^^and thus the convective blueshift is 
a significant effect and must be included in our model. 
We found that adding the convective blueshift changes 
the measurement of A by ^1 a, pushing it towards zero 
when the Vcb is included. 

We used a numerical model based on the work by Sh- 
porer & Brown (2011) similar to the approach used by 
(Albrecht et al. 2012). We made an initial assumption 
that the convective blueshift is similar to that of the sun 
to create a model grid for a range of Rp/ Ri, and impact 
parameters. We interpolated this grid at each step in 
the DE-MCMC chains. We left the velocity of the pho- 
tosphere (vcb) as a free parameter to account for the 
differences between HAT-P-17 and the sun. By defini- 
tion we expect vcb to be negative and for this reason 
we rejected models with positive vcb in the DE-MCMC 
chains. 



The semi-analytical model of Hirano et al. (2011) also 
includes three parameters that describe the sources of 
line broadening (/?, j^, and C)- Together with the rota- 
tional broadening of the star, these parameters provide 
a realistic analytical description of the observed line pro- 
files in the spectra. I3 includes both the Gaussian in- 
strumental line profile and the Gaussian dispersion from 
micro-tubulence. We adopted a fixed value of 3.0 km 
s~^ for /3 that is the result of summing in quadrature 
the width of the HIRES PSF (2.2 km s^^) and 2.0 km 
micro-turbulence broadening profile (Albrecht et al. 
2012). 'jH is the Lorentzian dispersion of the spectral 
lines primarily due to pressure broadening. We adopted 
a fixed value of 0.9 km s~^ that was found to be a good 
match to the HIRES spectra of several stars (Hirano et al. 
2011). The most significant of the line profile parameters 
is the macroturbulence broadening (C). We used equa- 
tion 3 from Valenti & Fischer (2005), 



c 



3.98 



cff " 



5770 K 



650 K 



km s 



(3) 



and Tcff — 5246 K (Howard et al. 2012) to calculate a 
value of 4.8 km s~^and assigned a conservative Gaussian 
prior of 3.0 km s"^ in the DE-MCMC analysis. We found 
that changing the prior centers on /3, 7, and C had little 
effect on the resulting posterior distributions of A, and 
for this reason we also could not remove the Gaussian 
prior on ^ without the DE-MCMC chains wandering into 
unphysical regions of parameter space. 

3.2.2. Convective blueshift 
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TABLE 3 

Rossitbr-McLaughlin MCMC results 



Parameter 



Value 



Units 



RM Step Parameters: 
log(P) 

Ttra 

i/e sin ui 
log(i?p/R.) 
log(a/i?*) 



cos J 

Ml + U'z 

V V sin i cos A 
V u sin i sin A 
/3 

7ir 
C 

VCB 

VCM 

S 

log (jitter) 



1.0144585 ±3.8e- 07 
2456165.8551 ±0.0011 

-U.044 _o.0068 

-0.214 ±0.015 
-0.9073 ±0.0032 

1 +0 0082 

l.OOOl _0.0085 
0123 +0 002S 
U.UIZO _0.0032 

736 +0-097 



log(days) 
BJDtdb 



0.687 


Vkm s 


-1 


24 +''■13 
'^■■^^ -0.2 


Vkm s 


-1 


4 1 +2-6 


km s 


-1 


= 0.9 


km s 


-1 


= 4.8 


km s 


-1 


-0.65 ±0.23 


log(km s~ 




11.13 i°:f 


m s 


-1 


-17.5 ilJ 


m s~iday 


-1 


0-289 1°:°^^ 


log(m s~ 





RM Model Parameters: 

P 10.3385231 



+9.1e-06 



days 



Ttra 


2456165.8551 ±0.0011 


BJDxDB 


e 


342 +0-0046 




201.5 tie 


degrees 




n lO-i^S +0.00092 
U.izolo -0.00091 






22.55 




i 


89.3 


degrees 




r^yr: +0.048 




vsini 


n +0.12 
0.56 14 


km s~l 


A 


19 

-16 


degrees 


/9 


4 1 +2.6 
^•^ -2.3 


km s~l 




= 0.9 


km s~l 


C 


= 4.8 


km s~l 


VCB 


-0.65 ±0.23 


km s~l 


VCM 


11.13 l°:f 


m s~l 


S 


-17.5 tti 


m s~iday-i 


jitter 


1 qc +0.38 
i.ao _o.34 


m s~l 



3.2.3. Results 

The results of the RM modeling are presented in Table 
3. Figure 3 shows the spectroscopic transit data with 
the best-fitting model overplotted and the resulting pos- 
terior distributions of wsinz^ and A. We measure the 
sky-projcctcd angle between the orbital angular momen- 
tum vector and the stellar rotation axis to be A = 191l\q 
degrees. This indicates that planet b's orbit is misaligned 
with the stellar rotation at a confidence level of only 1.2 
a. Our value of vsmi^, = 0.54 ± 0.15 km s~i is slightly 
larger than the value reported in H12 (vsinz^ = 0.3 ±0.5 
km s"'^), but well within the 1 a uncertainty from the 
SME analysis. 

We experimented with fixing vsm^^,, vcb, and the 
transit parameters and saw no significant changes in 
the resulting posterior distribution of A. When we ne- 
glect the convectivc blucshift in our model we measure 
a much more significant (presumably artificial) misalign- 
ment with A = 37 ± 12 degrees. We also examined the 



diagnostics from the Dopplcr analysis of the two outliers 
on cither side of the mid-transit. Wc found no evidence 
of systematic errors, poor fits, or other reasons to doubt 
the integrity of these model outlier points. Removing 
them from our fit did not change the results other than 
decreasing the reduced x^- 

3.3. Additional test for misalignment 

We also used the method of Schlaufman (2010) to in- 
dependently look for evidence of misalignment. This ap- 
proach compares the measured v sin i^. to an empirical 
estimate of the expected value of u = 2'KRi, / P-^ot, where 
Prot is the rotation period of the star based on the mass- 
age-rotation relations established from observations of 
the Hyades and Prasepe clusters summarized by Irwin 
& Bouvier (2009) . If the sky-projected inclination of the 
stellar rotation (z*) is close to 90° we would expect the 
measured vsini-^ to closely correspond to v. Since we 
know that the orbit of planet b is viewed nearly edge on 
{i = 89.3 

-0.17 degrees), an observed usini* significantly 
different than v would suggest spin-orbit misalignment. 
We use 



P^(M^,rO =P*,o(MO 



650 Myr 



1/2 



(4) 



from Schlaufman (2010) to calculate the expected rota- 
tion period of HAT-P-17 at the age given by H12. In 
the above equation, P*.o(-^*) is the rotation period of 
the star as a function of mass at an age of 650 Myr (12 
days), and — 7.8 ± 3.3 Gyr is the current age of the 
star. Assuming that the uncertainty in the age of HAT- 
P-17 is the dominant source of uncertainty, we calculate 
P,,(O.857M0,7.8± 3.3Gyr) = 42^% days. For = 90° 
we calculate v = 1.0]']q'2 km s~i. 

We compared our measured f sini^ = 0.56!tQ |^4 from 
the RM analysis to v by calculating the rotation statistic 
(0) from Schlaufman (2010) defined as: 



6 = 



t;sinz*obs 



(5) 



where wsini^^Qbg is measured from the RM analysis, and 
(Tt, and Cobs are the uncertainties on v and v sin ii, respec- 
tively. The difference between v sin i^, and v may suggest 
that the stellar rotation is inclined with respect to our 
line of sight. However, the value of = 1.9 is below 
the threshold for misalignment as defined by Schlaufman 
(9 > 2.9). This threshold accounts for the scatter in the 
empirical mass-age-rotation calibration which makes our 
determination of the initial rotation period of HAT-P- 
17 uncertain. The Schlaufman method provides modest, 
independent evidence of spin-orbit misalignment that is 
consistent with our obliquity measurement from the RM 
analysis. However, these two low significance measure- 
ments (each less than 2 a) do not conclusively show that 
the system is misaligned. 

3.4. Adaptive optics imaging 

We carried out high resolution and high contrast imag- 
ing with adaptive optics to check for near-by compan- 
ions in the context of understanding the architecture of 
the HAT-P-17 planetary system. Such companions are 
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Fig. 4. — Top: Keck/NIRC2 adaptive optics image. Bottom: 
contrast achieved based on the final reduced AO image. Our 
difFraction-Umited observations rule out the presence of compan- 
ions with Am > 7 mags for separations beyond 0.7". 



important in understanding the orbital evolution of the 
system. We find no evidence for off-axis sources in the 
immediate vicinity of HAT-P-17. To estimate our sen- 
sitivity to faint companions, we calculated the average 
contrast level achieved as a function of angular separa- 
tion. Specifically, we compared the peak stellar intensity 
to the standard deviation (ct) in scattered light within a 
square box of width 3 FWHM, where FWHM is the PSF 
full- width at half- maximum (also the size of a speckle). 
The standard deviation is evaluated at numerous loca- 
tions and the results are azimuthally averaged to create 
a contrast radial profile. 

We converted the contrast curve into a minimum de- 
tectable mass as a function of projected orbital separa- 
tion (figure 5, left panel) by interpolating the models of 
Baraffe et al. (2002) at the age and distance of HAT-P- 
17 from the analysis of H12. Assuming a circular orbit 
and Mp <C M^,, an order-of-magnitude approximation 
for the maximum RV slope caused by a fourth body in 
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Fig. 5. — Top: We convert the measured contrast to a mass- 
sensitivity curve using the Baraffe et al. 2003 evolutionary models. 
With an age of 7.8 Gyr, we could have detected any stellar com- 
panions (A/ > 80Mj) associated with HAT-P-17 at projected sep- 
arations beyond ^ 60 AU. Bottom: Same as top converted into a 
predicted RV slope using Equation 6. The region below the line is 
allowed by the data. The 3 a constraint on 7 from the RV analysis 
is 19 m s^^yr~^ which is just below the lower y-axis limit. 

the system is given by (Winn et al. 2009): 



GAIr sin ir 



7 ■ 



(6) 



We used this approximation to find the RV slope that 
would be produced by a planet at the minimum de- 
tectable mass as a function of projected orbital sepa- 
ration (figure 5, right panel). The minimum detectable 
mass at large separations is ~ 80A/j (coincident with 
the hydrogen-burning limit), far larger than the range 
of masses that are allowed by our RV fit and does not 
provide a good constraint on the orbit of planet c or a 
fourth companion. However, the RV data could still al- 
low for a long-period companion in a nearly face-on orbit 
or one that is currently near apsis which would minimize 
the radial velocity slope. The AO data help us rule out 
these scenarios for stellar/brown dwarf companions out- 
side - 50 AU. 

4. DISCUSSION 

HAT-P-17 is a rare planetary system with a transiting 
hot Jupiter and a long-period companion (HAT-P-17c). 
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There is no evidence for other companions. We have 
shown that the orbit of planet c is poorly constrained 
with the current RV data. We will not be able to conclu- 
sively measure the orbital parameters until a significant 
portion of the orbit has been observed. We find no evi- 
dence to suggest the presence of a massive 4th body. We 
modeled the RM effect of planet b and measure a possi- 
ble misalignment of the projected plane of the orbit and 
the rotation axis of the host star. 
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"1 I r 
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relative tidal-dissipation timescales Tce Tra 

Fig. 6. — Top: measured projected obliquity as a function of the 
alignment timescale calibrated from binary studies (Albrecht et al. 
2012). Stars with temperatures higher then 6250 K are shown with 
red filled symbols. Blue open symbols show stars with tempera- 
tures lower then 6250 K. Stars for which measured effective tem- 
peratures include 6250 K in their l-tr interval are shown by split 
symbols. We are computing the relative tidal dissipation rates as 
a function of stellar type, planet-to-star mass ratio, and orbital 
distance, using the scaling relations presented by Albrecht et al. 
(2012). Note that both timescales have been divided by 5 X 10^. 



Our constraints on a long-term RV trend (in addition 
to the two planet model) and the lack of companions seen 
in the adaptive optics image do not reveal any evidence 
for a 4th massive body in the system. However, given 
that the period of planet c reported by H12 was underes- 
timated we do not want to over- interpret any of planet c's 
parameters. Instead we urge the community to continue 
observing this interesting system in the coming years. 
We will be able to assess our measured 16.8 year orbital 
period in ~5 years when the RVs start to decrease rapidly 
as planet c approaches periapsis, and we will be able to 
start ruling out the short end of our estimate {Pc ~ 10 
years) in ~3 years. If it can be shown that the system is 
coplanar (this would require a spectacular observational 
effort by searching for transits of planet c) and apsidially 
locked then HAT-P-17 will be of further interest because 
it will give us a rare opportunity to probe the interior 
structure of an exoplanet by measuring the tidal Love 
number and quality factor through dynamical modeling 



(Batygin et al. 2009; Mardling 2010). 

We measure the sky-projected angle between the stel- 
lar spin axis and orbital angular momentum of the in- 
ner planet (stellar obliquity) by modeling the RM effect 
in Keck/HIRES RV data. The RM analysis suggests a 
slight spin-orbit misalignment of planet b with ~1.2 a 
confidence (A — 19li[g degrees). The Schlaufman method 
provides additional evidence for spin-orbit misalignment, 
but due to the dependence on somewhat uncertain stel- 
lar evolution models and the unknown initial angular mo- 
mentum of IIAT-P-17 we believe that the result from our 
detailed RM modeling to be more robust. 

Winn et al. (2010a) noted an emerging trend in the 
orbital obliquities of hot Jupiter hosts cooler than 6250 
K being generally well-aligned, and hot Jupiters around 
stars hotter than 6250 K having a wide distribution of A. 
IIAT-P-17 is a cool star around which we would normally 
expect to find spin-orbit-aligned hot Jupiters, but due to 
the relatively wide orbit of HAT-P-17b we do not expect 
that the planet would have had time to align itself with 
the spin of the host star (Winn et al. 2010a; Albrecht 
et al. 2012) if it were perturbed into a misaligned or- 
bit. This is consistent with our findings of a marginally 
non-zero A, however a more significant spin-orbit mis- 
alignment would provide stronger evidence against the 
disk-migration formation scenario. Figure 6 compares 
the timescales for alignment of this system with the sys- 
tems studied in Albrecht et al. (2012). If misaligned, 
this is the first multi-planet system in which a spin-orbit 
misalignment has been measured. We also note that our 
measurement of A is only slightly less likely to be consis- 
tent with zero, and coplanarity of planet c's orbit would 
provide strong evidence that the system migrated quies- 
cently in the protoplanetary disk. 

For a low-amplitude RM system like IIAT-P-17, we 
find that the convective blueshift is an important effect 
that must be included in our model for an accurate mea- 
surement of A. With 2-3 additional measurements of the 
RM effect we should be able to conclusively (^3 cr) de- 
termine if the system is misaligned which will help us 
understand the formation of the IIAT-P-17 system and 
other similar systems. 

We thank the many observers who contributed to the 
measurements reported here. We gratefully acknowledge 
the efforts and dedication of the Keck Observatory staff, 
especially Scott Dahm, Greg Doppman, Hien Tran, and 
Grant Hill for support of HIRES and Greg Wirth for 
support of remote observing. Finally, we extend special 
thanks to those of Hawai'ian ancestry on whose sacred 
mountain of Mauna Kea we are privileged to be guests. 
Without their generous hospitality, the Keck observa- 
tions presented herein would not have been possible. 
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